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1 Introduction 



Since it was shown that the path configuration of the random vicious walkers [1] is related with 
the Young tableaux [2-4], much attention has been paid on the statistical combinatorial prob- 
lems which are intimately related with the Young tableaux. Among them is the random permuta- 
tion [5], the random word [6], the point process [7, 8], the random growth model (the polynuclear 
growth model, oriented digital boiling model) [9, 10], the queuing theory [11], and so on. Inter- 
esting is that the scaling limits of these models have the universality that the fluctuation is of 
order A'''^'' with the mean being of order N, and that the asymptotic distribution of appropriately 
scaled variables is described by the Tracy-Widom distribution, which was originally identified 
with the limit distribution for the largest eigenvalue of the Gaussian unitary random matrix [12]. 
See Refs. 13-16 for a review. 

In this paper motivated from results in Ref. 17 and conjectures in Ref. 18, we introduce a 
physical model of the vicious walkers based on the hook Young tableaux. We shall study the 
scaling limit of certain probability, and clarify a relationship with the Tracy-Widom distribution. 

For our later convention we define the (M, A'^)-hook Schur functions (or, sometimes called the 
supersymmetric Schur function) [19] (see also Refs. 20-22), and denote some properties of the 
hook Young tableaux briefly. We set B = B+ U B_, and 

B+ = {6i, . . . , fM}, B_ = {Cm+I, • • • , fM+w}- (1-1) 

Hereafter we call i as positive (resp. negative) symbol when 6, e B+ (resp. e,- e B_). We fix an 
ordering in B as 

ei < e2< ■■■ < eM+N- (1-2) 

It should be noted that, though we use an ordering (1.2), following discussion can be applied for 
any other choices of ordering with |B+| = M and |B_| = A^. For a given Young diagram A, the 
semi-standard Young tableaux (SSYT) T is given by filling a number 1 , 2, . . . , M -i- in /I by the 
following rules; 

• the entries in each row are increasing, allowing the repetition of positive symbols, but not 
permitting the repetition of negative symbols, 
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• the entries in each column are increasing, allowing the repetition of negative symbols, but 
not permitting the repetition of positive symbols. 

We define the weight for SSYT T as 

wt(r) = m„6„, (1.3) 



where is the number of a's in T . Then the hook Schur function S x{.x,y) is given by 

SSYT T of shape A 

Here we have used 

Xi = e*"', for 6,- 6 B+, 

y.^^eM.j^ for 6m+; 6 B_. 

The Schur function s,i{x) in usual sense corresponds to a case of B_ = (D (N = 0), and the hook 
Schur function 5^ with B+ = (M = 0) reduces to the Schur function for the conjugate partition 
A'; 

s,(x,0) = sAx), SMy) = SA'(y). (1.5) 

Due to the rule of filling a number, the Young diagram A should be contained in the (M, A'^)-hook 
(see Fig. 1), and we have 

SA{x,y) = ^ Si,(x) SA'/^'iy). 
Furthermore when A contains the partition (A'^^), we have 

M N 

SAix,y) = Sf,(x) Syiy) ]~~[ ]~~[(x,- + yj), 
i=i j=i 

where the partitions fi and v are defined from A by yu, = Ai - N and vj = A'- - M, respectively. 
The Jacobi-Trudi formula helps us to write the hook Schur function as 

5,(^,3;) = det(c,,,y_,)^^..^^^^^^, (1.6) 
Here i(A) is the length of A, and c„ is given by the generating function, 

CO 

Hit;x)Eit;y) = YjC„f, (1.7) 

11=0 
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with 



E{t-y) = \\{\ + tyj). 
j 



(1.8) 



M 



^ N s- 

Figure 1 : (M, A^)-hook Young diagram must be contained in above "hook" region. 

This paper is organized as follows. In section 2 we introduce a model of vicious walkers as 
a generalization of the original model [1]. As far as we know, this model is first presented in 
this paper. We define the bijection from path configurations of vicious walker to the hook Young 
diagram. Especially we show a relationship between the length of the Young diagram and the 
number of movements of the first walker. This type of bijection was first given in Refs. 2, 3 for 
the original vicious walker model. In section 3 we give the probability of £{A) < £ m terms 
of the Toeplitz determinant. We further study the scaling limit of this probability based on the 
transformation identity from the Toeplitz determinant to the Fredholm determinant [23-25] in 
section 4. We apply the saddle point method to the Fredholm determinant following Refs. 10,26, 
and show that the scaling limit coincides with the Tracy-Widom distribution for the GUE [12]. 
In section 5 we consider some simple examples as a reduction of our model. Both the Meixner 
and the Krawtchouk ensembles can be regarded as a reduction of our vicious walker model. The 
last section is for conclusion and discussions. We briefly comment on the random word related 
with the hook Young tableaux. 

2 Vicious Walker 

We define a model of the random walkers which is related with the hook Schur function (1.4). 
The model is a generalization of one introduced in Ref. 1, and as will be clarified later an alge- 
braic property of the partition function is nothing but an identity in Ref. 17. 
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Evolution rule of vicious walkers is defined as follows. Initially there are infinitely many 
walkers at {...,-2,-1}, and we call each walker Pj whose initial point is -j. A walker is 
movable rightward if its right site is vacant. Walkers Pj+i,Pj+2, ■ ■ ■ are called successors of a 
walker Pj if they are next to each other in the order of the indices. We consider two types of 
time evolution (we assume that there are totally M + N time steps); first M-steps are referred as 
"normal" time evolution, and following A^-steps are as "super" time evolution. At a "normal" 
time evolution, a movable walker either stays or moves to its right together with an arbitrary 
number of its successors. Thus we draw 



move : 




On the other hand, at a "super" time evolution, a walker can move to its right any number of 
lattice units, though Pj cannot over-pass a position of Pj i at previous time. To realize this rule 
and to draw a non-intersecting path, it is convenient to depict this step as follows; 



arbitrary number of lattice 

Each path of the vicious walkers is required not to intersect each other. We see that the original 
model of vicious walkers [1] corresponds to a case of N = 0. After M + N-time steps, we 
denote Lj(n) as the number of right moves made by the walker Pj. Here n is a total number of 
movements of walkers. In Fig. 2 we give an example of path configuration of vicious walkers. 
In this case we consider totally 5-time steps (M = 3 and A^^ = 2), and the total step of right 
movements is n = 12 with (Li, L2, L3, L4) = (5, 4, 2, 1). 

It is now well known for the model of the original vicious walks [1] that we have the bijection 
from the path configuration of vicious walkers to the Young diagram [2]. This bijection can be 
easily generalized to our model as follows. For a path configuration (see, e.g.. Fig. 2), we draw 
Young tableaux A \- n with A'j = Lj(n). We insert in the 7-th column from top the times at which 
the j-th particle made a movement to its right. Notice that, for a super time evolution, we prepare 
number of time as many as lattice units the walker moved. For instance, in a case of Fig. 2, we 
put "1 2 4 4 5" in the first column, as Pi moves 2 lattice units rightward at time-4. Thus the top 
row is the list of times at which the walkers made their first movement, the second row is the list 
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Figure 2: Example of path configuration. For f < 3, a process is a "normal" evolution, while for f > 4 it becomes a 
super time evolution. 

of times at which the walkers made their second movement, and so on. It is clear that the normal 
time corresponds to the positive symbol B+ while the super time denotes the negative symbol 
B_ in SSYT. The evolution rule supports a consistency with ordering (1.2) in B, and we know 
that the map is indeed the bijection. Following this mapping the path configuration in Fig. 2 is 
mapped to SSYT given in Fig. 3. 
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Figure 3: Semi-standard (hook) Young tableaux with entries from B+ = {1, 2, 3) and B_ - {4, 5). 

To summarize, we have a one-to-one correspondence between the path configuration and 
SSYT; when n is the number of total moves of vicious walkers, we have A \- n, and the number 
Lj(n) of right movements made by the j-th walker is equal to the number of boxes in the j-th 
column of SSYT. Especially the length {(A) of partition coincides with Li(n). 



3 Partition Function and Toeplitz Determinant 



In the following we consider a model where, after totally n-step movements of the right movers, 
every walkers return to their initial positions by totally n-step left movements [3]. Here the 
number of normal (resp. super) time evolution is supposed to be Mi (resp. A'^i) in the first right 
moves, while the number of normal (resp. super) time evolution is M2 (resp. A^2) in the next 
left moves returning to their initial positions. The definition of normal and super time evolution 
in the left movers simply follows from that of right movers as a mirror image. Applying the 
bijection in previous section, the path configuration is denoted by pairs of SSYT A \- n, one is 
(Ml, A^i)-hook Young tableaux and another is (M2, A^2)-hook tableaux. 

We denote d,i(M, N) as the number of semi-standard Young tableaux of shape A with entries 
from B+ U B. (with |B+| = M and |B_| = A^). By definition, we have Si(?, . . . ,t,t, . . . ,t) = 

M N 

d^iM, N) t" for A \- n, and once the Young diagram A is fixed the number of SSYT d^iM, N) 
corresponds to the number of path configuration with fixed end points of right-moves. 

We have interests in the probability that the number of right movements of the first walker Pi 
is less than £, 



Here the probability "Prob" is defined as follows; we assign the weight t (we set < ? < 1) 
for every right- and left-moves, and regard the weight of totally n-step walk as t". Then each 
configuration of random walk, in which every walkers return to their initial positions after total 
2n-step, is realized with a probability f"/Z. An explicit form of the normalization factor Z 
will be given later. Based on the bijection map studied in previous section, we find that the 
probability (3.1) is given explicitly by 



{(A)<f 

Note that a normalization factor is set to be lim^^c^ Prob(Li <£) = !. 

To relate this probability with the random matrix theory, we follow a method in Ref. 6. 



Prob(Li < £). 



(3.1) 




(3.2) 
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Applying the Gessel formula to eq. (1.6), we have 



t{A)<C 

= 777^ fde n |e'^^-e'^f n//(e'^-x)£(e'^-3;)//(e-'^-z)£(e-'^-vv) 
= Dt{ip), (3.3) 
where ip{z) is defined by 

We have used Dc{(f) as the Toeplitz determinant for the function cp{z)\ Dc((p) is the determinant of 
£x£ matrix where an (/, j)-element is given by (fi^j with ip{z) = 2nez fn z". We note that eq. (3.3) 
was also given in Ref. 17. Thus our model of random walkers corresponds to a point process in 
Ref. 17 which was introduced as a generalization of Ref. 7. We note that the strong Szego limit 
theorem gives a generalization of the Cauchy formula, 

limDf(^)= — (3.5) 

i j,m n (1 " yj (1 ~ ^' 

We now apply a principal specialization ps which set jc, = aq' and yj = bqj [22]. In general 
we have 



ps(5'i(x,j)) = SA( aq,aq^,-- ; bq,bq^,... ) = q^-='''^- - 



and for a case of /} h n and (M, A'^)-hook Young diagram, by definition we have by setting 
a = b = t and q = I 

ps^^^^f.,^^i(Si(xi, ...,XM,yi,-- .,yN)) = dA(M,N) f. 
As a result, from eq. (3.3) we obtain the partition function as 

2 2 dAMuN,)dAM2,N2)t^" = D,icp), (3.6) 

n t(A)<i 
Ahn 



where 

(l+?z-i)^' {l+tzf ' 
(1 - ?z-i)^' (1 - tz)^'- 
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Due to the strong Szego limit theorem, we obtain a normalization factor Z as 



Z = limDK^) = (3.8) 



Combining these results, we get 



Prob(Li <^) = iz)K^). (3.9) 



4 Scaling Limit 

We study the asymptotic behavior of the probability (3.1). We note that in Ref. 18 the property 
of the scaling limit was conjectured. For our purpose, it is generally useful to rewrite the Toeplitz 
determinant with the Fredholm determinant. In fact, once we know the Toeplitz determinant, it 
is possible to rewrite it in terms of the Fredholm determinant [23-25]. Namely we have 

D,(^) = Zdet(l-DCa (4.1) 

where Z is defined in eq. (3.8), and %c is the matrix defined by 

DO 

Here a subscript denotes the Fourier component of the function, and we have used the Wiener- 
Hopf factor of(p,<fi = <fi+ (p^. 



(1-?Z)^2' (l_^^-l)M, 

Note that we have set < ? < 1. The probability (3.1) is now written by the Fredholm determi- 
nant as 

Prob(Li <£) = det(l - X(). (4.3) 

Using a representation (4.3) in terms of the Fredholm determinant, we study an asymp- 
totic behavior by applying the saddle point method following Refs. 10,26. We consider a limit 
Ma, Na ^ oo for a = 1 , 2 with fixed values; 

Ml A^i M2 

7r=mi, —=ni, — = m2. 

N2 N2 N2 
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In the Fredholm determinant (4.2), matrix elements are computed as 



2n\{l-tz)"^ {l+tiz) 

A path of integration in the former integral is chosen in a way that it surrounds z = -t, and that 
z = 1 /f is outside. On the other hand, a path of the latter integral includes both z = and z = t 
while it excludes z = -l/t. We set 

€ = cN2 + sn}, (4.4) 
where c is to be fixed later. For brevity, we define the function cr{z) by 

cr(z) = mi log(r -z) - Hi log(f + z) + log(l +tz)- m2log(l -tz) + (-mi + ni + c) logz. (4.5) 
Then above integrals are given by 

{^.l^-)-M-K-^ = (-1)^' ^lle^-(^)z^--^-^-^^'" ^ (-1)^' /i, 

(^-/^.)..,-.A..i =(-1)^^' ^lle-^-fe)z--^^'"-'-^-^^(-l)^'/2. 

We scale matrix indices as ii,j,k) {N^^^ x,N^^^ y,N^^^ w), and we consider to apply the 
saddle point method to integrals. 



r dz 

J 27ri 



^2 - 

e- 

in a limit A^2 oo- In these integrals, we fix a parameter c in eq. (4.4) so that we have a double 
saddle point, namely as a solution of a set of equations, 

do-iz) dV(z) 



dz dz^ 

i.e.. 



0, 



mi 1 m2 ni 

+ = c - m2 + 1 + h -, (4.6a) 



1-z/t 1+tz 1-tz l+z/t' 

nil Hi c - mi + Hi 1712 

j7^~(iTzf^ ? ^ (i-tzr~ (i + tzf 
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This set of equations is rewritten as 



t t tzo tzo „ . 

■ mi n\ - m2 - , (4.7a) 



where zo satisfies 



? - zo ?^ + zo 1 - ? Zo 1 + ? Zo ' 



+^ = ,^_L^ + ,^^. (4.7b) 



a-zo)2 a+zo)2 (i + ?zo)2 ii-tz^f 

We see that eq. (4.7b) always has a real solution in (-1 /?,-?) as far as Hi + because l.h.s. - r.h.s. 
of eq. (4.7b) changes from -oo to oo in z e (-1/f, -0- Generally real solutions of eq. (4.7b) are 
not only in (-1/?, -0, but to deform paths of integrals adequately we see that zo 6 (-1/^ -0 is 
a unique candidate of a double saddle point. For example, in a case of m\ = m2 and ni = I, real 
solutions of eq. (4.7b) are only z = ±1, and we can conclude that a double saddle point should 
be Zo = -1 from a discussion below. In a case of mi = > 1 and m2 = 1, real solutions of 
eq. (4.7b) are in {-\/t,-t), {t, l/t), (-oo,-l/r), and (1/?, oo) (the last 2 solutions exist only if 
mi = ni > l/f), and from a discussion to deform contours we see that only zo e (-1/^ -0 is 
possible as a double saddle point. Based on these cases, it may be natural to conclude that we 
choose Zo e (-!/?, -t) as a double saddle point. 

Hereafter we set a double saddle point zo so that zo 6 (-l/t, -t), and fix a parameter c by 
eq. (4.7a). With zo e {-l/t, -t), we find that c > from a definition (4.7a). With this choice of 
parameters, the fourth order equation (4.6a) has a real solution zo of multiplicity two, and other 2 
solutions are in {-t, t) and (1 It, oo). Around zo, we have a steepest descend path as in Fig. 4. As zo 
is a double saddle point, paths come into zo with angles +7r/3 and +2 7r/3. Following Ref. 10, we 
denote such paths as and Q~ respectively. We see that original paths explained above eq. (4.4) 
can be deformed smoothly to contours avoiding their singularities. Furthermore, we have 

1 dV(z) 



dz3 ^ _ 

Z-iO 

t /t-lzo t + 2zo l-2tzo l+2tzo 

^ \(F^ ^ (TTI^ " (1 - tzo? " (1 + tzo)' 

Zo ( l-t^ l+t^ It 
nil 7 Hi + 



l-tzo\{t-zof {t + zof (l + ^zo)^ 
where in the first equality we have used eq. (4.7a) to delete a parameter c, and in the second 
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equality we have erased m2 using eq. (4.7b). Recalling zo 6 (-l/t, -t) and Q <t < 1, we see that 

d^ a{z) 



dz3 



>0, 



(4.8) 



which shows that functions |e*^2 0-(zo)| j^^^g ^ maximum value at z = zo on a contour C*. 



0.5 



-0.5 




-2 







3 



Figure 4: Typical example of the steepest descent path 6^ is depicted as a bold line. Here we have set mi — 
«i - m2 - 1, and t - 1/2. A double saddle point is zo - -1, and other (simple) saddle points are 
(25 + 3 V4r)/16. We see that paths and come to a double saddle point zo at angles +n/3 and +2n/3 
respectively, and that another contour comes into (simple) saddle points with angle +n/2. A thin line 
denotes a local structure of the real part of the integrand around saddle points. 



With these settings, we have from the integral /i that 



^1/3 I ^N2 0-(z) N,^'^(w+y+s) 



r dz 

J 2;ri 

= A^y^ e^2 ""^""^ r e^ 2^2^^ (w+v+i) 
J 27ri 

J 2i;r' 

e+ 



'2 



cr cr 
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Here Ai(x) is the Airy function, 



Ai(z)= riie'(^-^/^\ 
J-oo 2n 



and we have set a parameter cr as 



o- = -Zo 





( 1 criz) 


\ 




[2 dz' 


z=zo^ 



1/3 



(4.9) 



We have cr > from eq. (4.8). In the same way, we have from the integral I2 that 



1/3 f _d£ 
J 27ri 



cr 



We then see that the kernel of the Fredholm determinant (4.2) is given by the Airy kernel, 



— dwAi( )Aii ) = - 

cr^ Jq cr cr cr Jn 



dwA2( + w)Ai( + w). 

cr cr 



As a result, we obtain 



lim Prob 



o-N. 



1/3 



< s 



(4.10) 



Here ^2(5) is the Tracy-Widom distribution [12] for the scaling limit of the largest eigenvalue of 
the Gaussian unitary ensemble, and is defined by 



F2(S) = det(l - ^Airy) 

dx (x - s) q{x) . 



(4.11) 
(4.12) 



The second equality is from Ref. 12, and q(x) is a solution of the Painleve II equation, 

q" = sq + 2q\ (4.13) 

with q(s) Ai(s) in ^ 00. 

A proof of convergence would be done along a line of Refs. 10, 26. 

To close this section, we comment on a case of ni =0. In this case, we further suppose 
that mi/f > 1 + m2. With this assumption, we see that there exists a real solution of eq. (4.7b) 
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in (-l/t,0). By setting this real solution to be zo 6 (-l/t,0), we can prove from eqs. (4.7a) 
and (4.9) that c > and cr > 0. Note that with this setting of a parameter c eq. (4.6a) has a real 
solution Zo of multiplicity two, and another solution is in {\/t, oo). See Fig. 5 as an example of a 
steepest descend path. As a result, we obtain the Tracy-Widom distribution (4.10) as a scaling 
limit. 




-2-10 1 2 3 



Figure 5: Example of the steepest descent path C'' for a case of «i = is depicted. Here we have set nii - 4, 1112 = 1, 
and t - 1/2. A double saddle point is zo = -0.68254, and there is a (simple) saddle point at 2.67684. As 
in Fig. 4, paths 6^ and 6" come to zo at angles +n/3 and ±2 n/3 respectively. Another contour comes into 
a (simple) saddle point with angle +7r/2, and it ends at t. A thin Une denotes a local structure of the real 
part of the integrand around saddle points. 

5 Some Special Cases 
5.1 Meixner Ensemble 

We consider a case, 

A/i = M2 = 0, i.e., nil = m2 = 0. 
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From the viewpoint of the random walkers, the vicious walkers move obeying only super time 
evolution rule. In this case the Toeplitz determinant (3.6) reduces to 

D,((l + -f ' (1 + tzf') = Z dA'mdAN2)t^" 



n C(A)<C 



J]J]d,(N,)d,(N2)f\ 



n fii<C 

where d^iN) denotes the number of (usual) semi-standard Young tableaux, and we have d^iN) = 
dxiN, 0) = d,\'{0, N) from eq. (1.5). The right hand side appeared in Ref. 7, and it gives an exam- 
ple of the discrete orthogonal polynomial ensemble as follows. Using the hook formula [22], 

\<i<i<M 

the r.h.s. gives 



d,{M)= W 



J-i 



(5.1) 



EE n 

n IJ[<( \l<i<j<N2 



j-i 



•l2 



n n 



V (=1 



j=N2'l 



j-i 



where we have assumed A^i > A^2- Introducing 

= l^j + N2- j, 

we obtain 

^j^' (M-A^2)! 
L y=o 

E 



(5.2) 



Prob(Li <£) = (l-tT' -t 



2\Ni N2 .~N2(N2-1) 



X 



jl(Ni-N2 + jV. 



ma.xlhi]<f+N2~i 



\<i<i<N2 



N2 



1=1 



N, -N2+h, 
hi 



f'\ (5.3) 



which is called the Meixner ensemble. 



In fact using the Borodin-Okounkov identity (4.2), the kernel of the Fredholm determinant 
can be written in terms of the Meixner polynomial 

(/ - j) j) = t'^^-^'-J (1 - r^)^'-^-^ ^ ^1 ■ (^') ■ (-A^2) 



\N2+j) \N2) 

x(Mj,^{i + N2\N,-N2 + \,t')-M„^__,{j + N2\N,-N2 + l,h 

- MM2-iii + N2; Ni-N2 + l,t^)- MN2U + N2; N1-N2 + 1, ?')), (5.4) 
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which has a well known form of the correlation functions of the random matrix (see, e.g., 
Ref. 27). Note that the Meixner polynomial is defined by 



b 



c 



A computation of the scaling limit can be done by the method in Section 4. A double saddle 
point, zo 6 -t), is explicitly solved as 

t+ ^Jn{ 

Zq = -- — - — =, 
1 + t yjnl 

and we obtain the Tracy-Widom distribution (4.10) with parameters c and cr defined by 

f(2V^ + (ni + 1)0 

C = —2 ' ^^-^^ 

cr = . (5.6) 

This result was derived by using the asymptotics of the Meixner polynomial in Ref. 7 (see also 
Ref. 28). 

5.2 Krawtchouk Ensemble 

We set 

Ni= M2 = 0, i.e., ni=m2 = 0, (5.7) 

The vicious walkers obey a rule of normal time evolution in the right moving, while they obey a 
rule of super time evolution in the left moving. In this case, eq. (3.6) is read as 



(1 + f7)^2 

^ z' ' n C{A)<{ 



2n 



This becomes the Krawtchouk ensemble [29] as follows (this type of the Toeplitz determinant 
was also studied in Ref. 10). When we substitute the hook formula (5.1) into above expression. 
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we see that the r.h.s. reduces to 



n 

7=0 



(Mi + j)! 



J 



X 



\}(Mj+N2-j)liM,+j-l-fij)\ 



By use of 



f^j =l^j + N2- j. 



this gives 

Prob(Li <i) = (l+ h'''' 



U 7!(A^2 + Mi-1)! 



X 



E 



max|/!,j<^+A'2-l 



\<i<j<N2 



(5.8) 



which is the Krawtchouk ensemble. 



The kernel of the Fredholm determinant is computed explicitly from eq. (4.2), and it is given 
in terms of the Krawtchouk polynomial as a form of the correlation functions; 

(i-])%{i,i) = . (Ml + A^2 - 1) \^ \ ■ \, . 



X Lfi:w,(/ + A^2; -; Ml + A^2 - 1) • KN2-1U + N2; 



f- 



.Ml +A^2 - 1) 



2 2 

- i^^,_i(/ + A^2; 7^, Ml + A^2 - 1) • is:A^2(7 + A^2; M,+N2- 1)). (5.9) 

Here the Krawtchouk polynomial is defined by 



K„(x;p,N) = 2Fy 



-n, -X 
-N ' 



We note that we have 



Kn{x-p,N) = M„{x; -N,-^). 

P- 1 



The scaling limit is also computed by the saddle point method [10]. In this case we suppose 

nil > t^, and we have a double saddle point zo ^ {-I It, 0) as 

ymj" + 1 



zo = 
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and we obtain the Tracy-Widom distribution (4.10) where parameters c and cr are defined from 
eqs. (4.7a) and (4.9) as 

f(2V^ + (mi - 1)0 
c= — , (5.10) 

ri/3(VJi7-o2/3(i+rV^)2/3 
a = — . (5.11) 

One sees that this result coincides with that of Ref. 29 derived by use of asymptotics of the 
Krawtchouk polynomial. 

5.3 Symmetric Case 

We consider a case, 

Ml = M2, Ni = N2, (i.e., nil = m2 = m, rii = 1), (5.12) 

namely in right- and left-movements we have equal number of normal and super time evolu- 
tions. Unfortunately we are not sure whether this model is related with the discrete orthogonal 
ensemble, but the parameters of the scaling function can be simply solved as follows. 

In a scaling limit N2 00, we obtain the Tracy-Widom distribution (4.10) by applying the 
saddle point method. In this case a double saddle point is zo = -1> and parameters in eq. (4.10) 
are computed from eqs. (4.7a) and (4.9) as 

2?(1 +m + (l -m)0 

c = Y~p. ' ^^-^^^ 

^= —2 • (5-14) 

6 Conclusion and Discussion 

We have introduced a generalization of the vicious walker model in Ref. 1 . We find that there 
exists a bijection map between the path configuration of vicious walkers and the hook Young 
diagram as in the case of the original vicious walkers. We have exactly computed a probability 
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that the number of right movements of the first walker is less than i, and have given a formula in 
terms of the Toeplitz determinant. We have further studied a scaling limit of the probability based 
on the Borodin-Okounkov identity which relates the Toeplitz determinant with the Fredholm 
determinant, and have obtained the Tracy-Widom distribution for the largest eigenvalue of the 
Gaussian unitary random matrix. Other models which belong to the orthogonal or the symplectic 
universality classes are for future studies. 

In the case of the vicious walker model, crucial point is that there exists the bijection map 
from the path configuration to a pair of the semi-standard (hook) Young tableaux. As was well 
studied [6], a pair of SSYT and the standard tableaux is related with the problem of the random 
word. We can define the model of the random word which is related with the hook Young diagram 
as follows [30] . We consider a random word by choosing from a set B+ uB_ with B+ = { 1 , . . . , M} 
and B_ = {M + I, . . . ,M + N}. When a word of length n is given, we have a generalization of 
the Robinson-Schensted-Knuth (RSK) correspondence [31,32] (see also Ref. 33 for invariance 
under ordering of symbols); we have a bijection between a word of length n and pairs (P, Q) of 
tableaux of the same shape A\- n{Ph SSYT from B, and the recording tableaux Q is the standard 
Young tableaux). Rule to construct pairs of tableaux is essentially same with the original RSK 
correspondence (see, e.g., Ref. 21,22), and a difference is only that negative symbols can bump 
himself while positive symbols cannot bump himself. Then for random word with length n, the 
probability that the length of longest decreasing (strictly decreasing for positive symbols while 
weakly decreasing for negative symbols) subsequence is less than or equal to I is then given by 



where is the number of standard Young tableaux. 

This can be rewritten in terms of the Toeplitz determinant based on eq. (3.3). We use the 
exponential specialization [22], 




(6.1) 



e(A)<e 



(6.2) 



where the power sum symmetric function p„ is given by 
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Acting on the hook Schur function, we have 



tx{S,(x,y))=f-, 



for A \- n. By applying the exponential specialization to (x,y) and the principal specialization 
ps„^^^,.^^i to (z, w) in eq. (3.3), we get 



As a consequence, the Poisson generating function of the probability (6.1) is given by the Toeplitz 
determinant of function O. As was seen from the fact that the kernel (6.4) can be given from 
eq. (3.7) as an appropriate limit, the scaling limit of eq. (6.3) reduces to the Tracy-Widom distri- 
bution as was shown in Ref. 16 for a case of = 0. Detail will be discussed elsewhere. 

It was shown in Ref. 6 that the generating function (6.3) with N = have an integral repre- 
sentation in terms of solutions of Painleve V equation. It remains for future studies to clarify a 
relationship between the Toeplitz determinant (6.3) in a case ofNi^O and the Painleve equations, 
especially integral solutions of the Painleve equation given in Ref. 34. 

Note Added: After submitting this paper, Ref. 35 appeared on net. Therein studied was a limit 
theorem of the "shifted Schur measure", where the probability is defined in terms of the Schur 
2-functions [20]. To apply an a method of Ref. 10 they obtained the Fredholm determinant 
after a finite perturbation of a product of Hankel operator, but their main result on a scaling limit 
exactly coincides with our results (4.10) with Mi = A'^i and M2 = N2 (subsequently one sees 
that their result for r = 1 coincides with our above results (5.13)- (5.14) with m = 1). This 
coincidence may originate from a property of the Schur Q-function. The Schur Q-function is 
defined by filling "marked" and "unmarked" positive integers to the shifted Young diagram; a 
rule of filling these numbers is much the same with a rule for the semi-standard hook Young 
tableaux explained in Introduction, once we identify unmarked (resp. marked) numbers with 
positive (resp. negative) symbols. It will be interesting to investigate this connection in detail. 




(6.3) 



where 



0(z) = e' 



(1-z) 



(6.4) 
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